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equations for the radial parts of the spinor wave function are shown to involve a squared 
Dirac operator for the free case, whose essential self-adjointness is proved by using the Weyl 
limit point-limit circle criterion, and a 'perturbation' resulting from the potential. One 



Qh! then finds that a potential of Coulomb type in the Dirac equation leads to a potential term 



in the above second-order equations which is not even infinitesimally form-bounded with 
respect to the free operator. Moreover, the conditions ensuring essential self-adjointness 
of the second-order operators in the interacting case are changed with respect to the free 
case, i.e. they are expressed by a majorization involving the parameter in the Coulomb 
potential and the angular momentum quantum number. The same methods are applied 
to the analysis of coupled eigenvalue equations when the anomalous magnetic moment of 
the electron is not neglected. 



1. Introduction 

In the same year when Dirac derived the relativistic wave equation for the electron [1], 
the work of Darwin and Gordon had aheady exactly solved such equation in a Coulomb 
potential in three spatial dimensions [2, 3]. Since those early days, several efforts have 
been produced in the literature to solve the Dirac equation with other forms of central 
potentials, until the recent theoretical attempts to describe quark confinement [4-8]. In 
the present paper we study the mathematical foundations of the eigenvalue problem for a 
massive spin-i field in a central potential V{r) on R^, without specifying a priori which 
function we choose for V{r). In other words, we prefer to draw conclusions on V{r) from 
a careful mathematical investigation. 

By doing so, we hope to elucidate the general framework of relativistic eigenvalue 
problems on the one hand, and to develop powerful tools to understand some key features 
of central potentials on the other. For this purpose, in section 2 we focus on the radial parts 
of the spinor wave function, casting the corresponding second-order differential operators 
in a convenient form for the subsequent analysis. In section 3, the Weyl limit point-limit 
circle criterion [9] is used to prove that the squared Dirac operator for the free problem is 
essentially self-adjoint on the set C^{0,oo) of smooth functions on (0, oo) with compact 
support away from the origin. In section 4 some boundedness criteria for perturbations [9, 
10] are first described and then applied when the potential in the original Dirac equation 
consists of terms of Coulomb and/or linear type. The effects of the anomalous magnetic 
moment of the electron are studied in section 5. Concluding remarks and open problems 
are presented in section 6. 

2. Second-order equations for stationary states 

For a charged particle with spin in a central field, the angular momentum operator and 
the parity operator with respect to the origin of the coordinate system commute with the 
Hamiltonian. Thus, states with definite energy, angular momentum and parity occur. The 



corresponding spinor wave function for stationary states reads [11, 12] 



^ = 






where ^j,i,m and flj^i'^m are the spinor harmonics defined, for example, in [11, 12], and 
l = j±lJ' = 2j-l. 

The stationary Dirac equation in a central potential V{r) takes the form (toq being 
the rest mass of the particle of linear momentum p) 



a • p —niQC^ + Vir)} \x 
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and leads eventually to the following coupled system of first-order differential equations 
(having defined F{r) = rf{r) and G{r) = rg{r)): 



-^ + -] G{r) = (Ai - W{r))F{r) 
dr r 



(2.1) 



CI rC 1 ._ , > 

_ + _ F{r) 
dr r 



{-X2-W{r))G{r) 



(2.2) 



where k = —I — 1 (if j = / + ^) or / (if j 



|), and we have defined 



W{r) 



V{r) 

he 
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Ai = 



A2 = 



E + rriQC^ 

he 

—E + rriQC^ 
he 
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Equation (2.1) yields a formula for F{r) which, upon insertion into Eq. (2.2), leads to the 
second-order equation 



-T^+p{r)— + q{r) 
ar^ dr 



G{r) = 



(2.6) 



where 

W'(r) 
(Ai -W{r)) 

q{r) = _M^±11 + tp^r) + W\r) + (A2 - Xi)W{r) - A1A2. (2.8) 

Equation (2.6) should be supplemented by the boundary condition G{0) = 0. It then 
describes a Sturm-Liouville equation non-linear in the spectral parameter. In [13], the 
equivalence has been proved of the radial Dirac equations (2.1) and (2.2) to the parameter- 
dependent Sturm-Liouville equation (2.6) (the parameter A used in [13] corresponds to our 
E, and moc^ = 1 units are used therein). By equivalence we mean that, under suitable 
assumptions on the potential, the function G solving Eq. (2.6) is found to belong to the 
prescribed space ifd(R_|_), i.e. the space of absolutely continuous functions on [0, 00) which 
are square-integrable on R_|_ jointly with their first derivative and vanish at the origin. Now 
we can use a well known technique to transform Eq. (2.6) into a second-order equation 
where the coefficient of -^ vanishes. This is achieved by defining the new function O such 
that [14] 

0(r) = G{r) exp - / p{r)dr. (2.9) 

In the few cases where exact analytic formulae are available in the literature one studies 
indeed Eq. (2.6) and its countepart for F (see Eq. (2.13)). However, Eq. (2.9) has 
the advantage of leading to a second-order equation for O in a form as close as possible 
to 'perturbations' of Schrodinger operators, and is hence preferred in our paper devoted 
to qualitative and structural properties. In non-relativistic quantum mechanics, such a 
method leads to a unitary map [9] transforming the radial Schrodinger equation in a 
central potential into an equation involving a radial Schrodinger operator —-j^ + U{r) 
acting on square-integrable functions on R_|_ which vanish at the origin. In our relativistic 
eigenvalue problem the transformation of the Hilbert space of square-integrable functions is 
no longer unitary, but remains of practical value. All non-linear properties of the resulting 
Sturm-Liouville boundary- value problem are in fact encoded into a single function playing 



the role of parameter-dependent potential term (see below), rather than two functions p 
and q as in (2.6)-(2.8). The function O is then found to obey the differential equation 



d^ l{l + l) 



dr"^ 
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(r) = -W'^{r) + 
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Such an equation may be viewed as follows: since the potential W 'perturbs' the 'free' 
problem for which W vanishes in Eqs. (2.1) and (2.2), in the corresponding second-order 
equation (2.10) one deals with a 'free operator' 



A 



^\'l±ll!or.n 



dr"^ 



r" 



0,1,... 



(2.12) 



perturbed by the multiplication operator Pw,E{r) defined in (2.11). An interesting pro- 
gramme is therefore emerging at this stage: 

(i) First, prove (essential) self-adjointness of the 'free' operator Aj, on a certain domain. 

(ii) Second, try to understand whether the operator A^^ + Pw,E{f) in Eq- (2.10) remains 
self-adjoint on the same domain. If this condition is too restrictive, try to derive all 
properties of this 'perturbed' second-order operator. 

If one first uses Eq. (2.2) to relate G{r) to ^ and F, one finds instead the Sturm- 
Liouville equation (cf [13]) 



^2+P(r)^^+^(r) 



F{r) = 



(2.13) 



supplemented by the boundary condition -F(O) = 0, where (cf (2.7) and (2.8)) 



p{r) 



W'{r) 



{\2 + W{r)) 



(2.14) 



q{r) = - ^^^^^ ^^ - ^p{r) + W\r) + (A2 - Xi)W{r) - A1A2. (2.15) 






Thus, after defining (cf (2.9)) 



0(r) = F(r) exp - / p{r)dr (2-16) 



one finds for 0(r) the second-order differential equation 



^2 k{k - 1 



+ + PwMr) 



2 ^2 ^4^1 



0(r) = ^^^,y- 0(r) (2.17) 



1 W" ?,fW'^'^ 



having now defined (cf (2.11)) 

Since k — —I — 1 \i j = I + ^, and k — I ii j = I — ^^ the 'free' operator in Eq. (2.17) reads 
now 



f_ {l + l){l + 2) 
dr'^ j,2 



yi; = --— + ^ '; ' for all / = 0, 1, ... (2.19a) 



^^2 j,2 



Al = --j-^+ ^ ^^ ' for all / = 1, 2, ... . (2.196) 



Note that Pw,e{t) has a second-order pole at Ai = VT (see (2.11)) and Pw,E{'r) has a 
second-order pole at A2 = —W. Thus, the analysis of the interacting case (i.e. with 
W{r) 7^ 0) is performed in section 4 at fixed values of E and away from such singular 
points. 

3. Weyl criterion for the squared Dirac operator in the free case 

The self-adjointness properties of the free operator (2.12) should be studied by considering 
separately the case / > and the case / = 0. For positive values of the quantum number 
/, Aj, turns out to be essentially self-adjoint. This means, by definition, that its closure 

6 



(i.e. the smallest closed extension) is self-adjoint, which implies that a unique self-adjoint 
extension of Al. exists [15]. In general, if several self-adjoint extensions exist, one has to 
understand which one should be chosen, since they are distinguished by the physics of the 
system being described [9, 15]. This is why it is so desirable to make sure that the operator 
under investigation is essentially self-adjoint. We here rely on a criterion due to Weyl, and 
the key steps are as follows [9]. 

The function V is in the limit circle case at zero if for some, and therefore all A, all 
solutions of the equation 



are square integrable at zero, i.e. for them 



(/?(x) = \(p{x) (3.1) 



\^{x)\^dx<oo (3.2) 







with finite values of a, e.g. a g]0, 1]. If V{x) is not in the limit circle case at zero, it is 
said to be in the limit point case at zero. The Weyl limit point-limit circle criterion states 
that, if F is a continuous real- valued function on (0, oo), then the operator 

0^-^,+V{x) (3.3) 

is essentially self-adjoint on C^(0, oo) if and only if V{x) is in the limit point case at both 
zero and infinity. The property of being in the limit point at zero relies on [9] 

Theorem 3.1 Let V be continuous and positive near zero. If 

V{x) > jx-^ (3.4) 

near zero, then O is in the limit point case at zero. 

The limit point property at oo means that the limit circle condition at oo is not 
fulfilled, i.e. the condition 

\ip{x)\'^dx < oo (3.5) 

7 



does not hold. To understand when this happens, one can use [9] 

Theorem 3.2 If V is differentiable on (0, oo) and bounded above by a parameter K on 

[1, oo), and if 

"^^ oo (3.6) 



'i ^K - V{x) 

V'{x)\V{x)\~^ is bounded near oo (3.7) 

then V{x) is in the hmit point case at oo. 

Thus, a necessary and sufficient condition for the existence of a unique self-adjoint exten- 
sion of O is that its eigenfunctions should fail to be square integrable at zero and at oo. 
Powerful operational criteria are provided by the check of (3.4), (3.6) and (3.7), which only 
involve the potential. 

In our problem, for all / > 1, the 'potential' Vi{r) = ^ 2 is of course in the limit 
point at zero, since the inequality (3.4) is then satisfied. Moreover, Vi{r) is differentiable 
on (0, 00), bounded above by xi — K^ + 1) ^^ [l^ '^^)^ ^^^ such that 

-MX = 00 (3.8) 



y^;r^ ^JimyA -/^^^i 



V,'{r)mr)n = - ,-,—- foi- all r. (3.9) 

Hence all conditions of theorem 3.2 are satisfied, which implies that Vi{r) is in the limit 
point at oo as well. By virtue of the Weyl limit point-limit circle criterion, the free operator 
Al defined in (2.12) is then essentially self-adjoint on 0^(0, oo) for all / > 0. 

When / = 0, however (for which k = —1), A^ reduces to the operator —-§^: which 
has deficiency indices (1, 1). Recall that for an (unbounded) operator B with adjoint B\ 
deficiency indices are the dimensions of the spaces of solutions of the equations B^u = ±iu. 
More precisely, one defines first the deficiency sub-spaces {D{B'^) being the domain of B^) 

n+{B) = {ueD{B^):B^u = iu} (3.10) 



n-{B) = {ueD{B^):B^u = -iu} (3.11) 

with corresponding deficiency indices 

n+{B)= dim n+{B) (3.12) 

n-{B) = dim n-{B). (3.13) 

The operator B is self-adjoint if and only if n+{B) = n-{B) = 0, but has self-adjoint 
extensions provided that n+{B) = n-{B) [9, 15]. In our case, half of the solutions of 

the equations (^4^) u = ±iu are square-integrable on R_|_, which implies that n+(74^) = 
n_(A^) = 1. This is easily proved because such equations with complex eigenvalues reduce 
to the ordinary differential equation [9] 

e"^ = ie"'^ (3.14) 



dr 
and 



d"^ 

e-- = -ie-^ (3.15) 



dr'^ 



In the former case, on setting a = pe'^, with p and 6' G R, one finds p = ±1,6 = — f , 
which leads to the two roots of the equation — a^ = i: 

1 i 

In the latter case, u solves the algebraic equation uj^ = i, and hence one finds the roots 



Only the roots a2 and uj2 are compatible with the request of square-integrable solutions 
of (3.14) and (3.15) on R_(_, and hence one finds nj^{A^) = n_(A^) = 1 as we anticipated. 
This property implies that a one-parameter family of self-adjoint extensions of A° exists, 
with domain D{A^) given by 

D(^°) = {ue L^{R+) :u,u' e ACioc(R+);w" e L^{n+); 
u{0) = f3u (0)} . (3.20) 

Here ACioc(R'+) denotes the set of locally absolutely continuous functions on the positive 
half-line, the prime denotes differentiation with respect to r, and /3 is a real-valued param- 
eter. Bearing in mind the limiting form of Eq. (2.10) when / = and W = 0, this means 
that one is studying the case characterized by 

,.(£!_^<0 (3.21) 

for which the square-integrable eigenfunction of —-;^ reads {a being a real constant to 
ensure reality of E) 

u{r) = a e-'^VW _ (3.22) 

On defining 



(u,v)= / u*{r)v{r)dr 
Jo 

the boundary condition in (3.20) is obtained after integrating twice by parts in the integral 

defining the scalar product (A'^u^v) to re-express it in the form f m, {A^)'^vy with u in 

the domain of A^ and v in the domain of the adjoint {A'^Y . One then finds that both u 
and V should obey the boundary condition (3.20). In the light of (3.20)-(3.22) one obtains 
the very useful formula 

1 = -(3^/\X\ (3.23) 



which implies 



E' = mlc'-^. (3.24) 
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This means that in a relativistic problem a lower limit for 0^ (and hence for |/3|) exists, to 
avoid having E^ <Q. 

To complete the analysis of squared Dirac operators in the free case, one has also to 
consider the operators Aj, defined in (2.19a) and (2.19b). The former has a 'potential' term 
- — 4 — which is in the limit point case at both zero and infinity for all/ > 0. The latter 
has a 'potential' term ^ ^ ' which is in the limit point at zero with the exception of the 



7 JZ 

value 1 of the quantum number /, for which A\. reduces to the operator — -t-^, and hence 



we repeat the logical steps proving that such an operator has a one-parameter family of 
self-adjoint extensions. Once more, their domain is given by Eq. (3.20). 

4. Second-order operators in the interacting case 

Now we would like to understand whether the general results on perturbations of self- 
adjoint operators make it possible to obtain a better understanding of effects produced by 
the central potential W{r) in Eqs. (2.10) and (2.17) (the essential self-adjointness of the 
Dirac Hamiltonian with non- vanishing W is studied in [16], and several comments can be 
found in the following sections). For this purpose, the key steps are as follows [9]. 

(i) Let A and B be densely defined linear operators on a Hilbert space H with domains 
D{A) and D{B), respectively. If D(A) C D(B) and if, for some a and b in R, 

\\Bip\\<a\\Aip\\+b\\ip\\ iov all ip e D{A) (4.1) 

then B is said to be A-bounded. The infimum of such a is called the relative bound of B with 
respect to A. If the relative bound vanishes, the operator B is said to be infinitesimally 
small with respect to A. 

(ii) The Kato-Rellich theorem states that if A is self-adjoint, B is symmetric, and B is 
A-bounded with relative bound a < 1, then A + 5 is self-adjoint on D{A). 

(iii) If the potential V can be written as 

V = Vi + V2 (4.2) 
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with Vi e L^(R^) and V2 G L°°(R^), and ifV is real-valued, then the operator —A+V{x) 
is essentiaUy self-adjoint on Co°(R"^) and self-adjoint on D(—A). As a corollary, the 
operator — A — ^ is essentially self-adjoint on Co°(R^). 



(iv) An analogue of the Kato-Rellich theorem exists which can be used to study the case 
when B is not A-bounded. The result can be stated after recalling the following definitions. 
Let A be a self-adjoint operator on H. On passing to a spectral representation of A 
with associated measures {/in}^=i on the spectrum of A, so that A is multiplication by x 
on the direct sum ©^^2^L^(R, dfXn), one can consider 






N 
^ = i {'^n{x)y^^l : "^ I \x\\lpn{x)\''d^n < 00} (4.3) 

and hence define, for ip and ip eT, 

li^,'^) = ^ XCp*^{x)^n{x)dfin- (4.4) 

Such a g is called the quadratic form associated with A, and one writes 

Q{A)=X. (4.5) 

The form domain of the operator A is then, by definition, Q{A), and can be viewed as the 
largest domain on which q can be defined. 

(v) The KLMN theorem states that, if A is a positive self-adjoint operator and if l3{(f, il)) 
is a symmetric quadratic form on Q{A) such that 

\(3{ip, ip) I < a((^, Aip) + h{ip, if) for all p e D{A) (4.6) 

for some a < 1 and 6 G R, then there exists a unique self-adjoint operator C with 

Q{C) = Q{A) (4.7) 
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and 

((^, Cij) = {ifi, Atjj) + (3{ip, ij) for all </?, ^ G g(C) . (4.8) 

Such a C is bounded below by —b. 

(vi) If ^4 is a positive self-adjoint operator, and S is a self-adjoint operator such that 

Q{A) c Q{B) (4.9) 

and 

\{if,By^)\<a{ifi,Ay^) + b{^,^) for all v? e D(^) (4.10) 

for some a > and 6 G R, then B is said to be relatively form-bounded with respect 
to A. Furthermore, if a can be chosen arbitrarily small, B is said to be infinitesimally 
form-bounded with respect to A. 

(vii) If the operator B is self-adjoint and relatively form-bounded, the parameter a being 
< 1, with respect to a positive self-adjoint operator A, then the KLMN theorem makes 
it possible to define the 'sum' A + S, although this mathematical construction may differ 
from the operator sum. In particular, B can be form-bounded with respect to A even 
though the intersection of their domains may be the empty set. 

(viii) The KLMN theorem is physically relevant because it leads to the definition of Hamil- 
tonians even when the Kato-Rellich criterion is not fulfilled. In other words, the request of 
dealing with L'^+L'^ potentials is too restrictive. For example, the potential Voi{r) = — r~" 
belongs to L^ -{- L°° only if a < |. However, if a G [f , 2), one can use the KLMN theo- 
rem because, for all a < 2, one can prove that — r~" is infinitesimally form-bounded with 
respect to —A [9]. 

In our problem, the 'potential' terms in Eqs. (2.10) and (2.17) are given by (2.11) and 
(2.18), respectively. If the potential W{r) is of Coulomb type, i.e. (7 being a dimensionful 
constant) 

W{r) = - (4.11) 
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the singular behaviour of Pw,E{r) as r ^ is dominated by (for a fixed value of E) 



il' + l + k) 



and the singular behaviour of Pw,E{'r) as r ^ is given instead by (again for a fixed value 

of ^) 

(7^ + i - k) 



Thus, as r ^ 0, the operators on the left-hand sides of both (2.10) and (2.17) reduce to 



Ijj- 






dr 



(4.12) 



In the operator L^, the coefficient of r~^ is no longer greater than or equal to | (see 
(3.4)) for the same values of / ensuring essential self-adjointness of the free problem. The 
inequality 



k' 



7 



1>5 

4-4 



is instead fulfilled by 



if /c = — / — 1, and by 



{l + lf >7^ + l forall/ = 0,l,... 



/^ > 7^ + 1 for all/ = 1,2, 



(4.13) 



(4.14) 



(4.15) 



a k = I. Our result implies that, for all \k\ > 2, essential self-adjointness on Cq°{0, oo) of 
the second-order operators on the left-hand sides of (2.10) and (2.17) is obtained provided 
that I7I < Vo. This reflects the fact that a Coulomb potential in the flrst-order system (2.1) 
and (2.2) leads to 'potential' terms in the second-order equations (2.10) and (2.17) which 
are not even inflnitesimally form-bounded with respect to the squared Dirac operators 
in the free case, because both the potential terms and the free operators contain terms 
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proportional to r ^. To study the limit point condition at infinity, we try to majorize the 
'potential' Pw,e obtained from the Coulomb potential (4.11), and we find that 



|Pv7,E(r)|<2^ + 



3„2 



|7(l + /^)|(Ai + |7|) + |7 



he xj 

if r G [1,00). Moreover, the integral (3.6) diverges when V is replaced by Pw,Ej and the 
condition (3.7) is fulfilled as well, because 

3 1 

The check of (3.6) and (3.7) for Pw,e leads to the same results, and hence we use the Weyl 
criterion of section 3 to conclude that, for fixed values of E, essential self-adjointness on 
Cq°{0, 00) of the second-order operators in Eqs. (2.10) and (2.17) holds provided that the 
inequality A;^ — 7^ > 1 is satisfied. This rules out / = in (4.14) and ^ = 1 in (4.15). One 
then finds that I7I < a/3 as we said before. 

The limiting form (4.12) is not affected by the addition of parts linear in r [6, 17, 18] 
to the right-hand side of (4.11), because the singular behaviour of Pw,E{r) at fixed values 
of i? as r ^ is still dominated by the Coulomb potential. By contrast, a purely linear 
potential 

W{r) = Tr (4.16) 

satisfies the request of infinitesimal form-boundedness of Pvk,e(^) with respect to the 
squared Dirac operators in the free case, because then the singular behaviour of Pw,E{r) 
as r ^ is expressed by — -^^ and the singular part of Pw,e{'>') as r -^ reads —Y~r- 
However, one might consider linear terms with compact support, i.e. vanishing for all r 
greater than some finite ro, or weighted with exponential functions which ensure a fall-off 
condition at infinity, e.g. the potential (cf [17]) 

W{r) = - + rre-'''' (4.17) 
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where fi is positive. In such a case, the hmiting behaviours of Pw,e as r ^ and as 
r ^ oo are stiU dominated by the Coulomb part in the potential W, and hence we find 
again essential self-adjointness on Cq°{0, oo) provided that /c^ — 7^ > 1. 

In the physical literature, however, the potential has not been written in the form 
(4.17). To achieve quark confinement, a purely linear term has instead been added to 
the Coulomb part, considering also a split of the additional part into Lorentz scalar- type 
and Lorentz vector-like potentials. Furthermore, such a vector contribution is sometimes 
omitted in a phenomenological analysis, bearing in mind its non-perturbative nature (since 
the perturbative part has instead vector nature) [19]. Needless to say, such arguments are 
not compelling. 

5. Inclusion of the anomalous magnetic moment 

The second-order operators that we have analyzed in the interacting case (see again Eqs. 
(2.10) and (2.17)) are not 'squared Dirac operators' because the eigenvalues of the Dirac 
operator occur in their 'potential term'. It is therefore important to compare more carefully 
the predictions of the second-order equation for O (and O) with the results obtained from 
squared Dirac operators studied in [16]. The latter are used in [16] because a theorem 
ensures that, given the (abstract) Dirac Hamiltonian 

if one of the operators D-D^ or Dj^D- is essentially self-adjoint, then the operator T is 
essentially self-adjoint as well, where W+ and W- take into account the rest mass and the 
potential (see theorem 5.9 in [16]). 

Let us now consider the effect of the anomalous magnetic moment /i of the electron in 
a central potential V{r). With the notation of our section 2, the resulting set of coupled 
eigenvalue equations is found to be [16] 



ar r 



G{r) = (Ai - W{r))F{r) (5.2) 
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-^ + --fiW'{r) 
dr r 



F{r) = {-\2-W{r))G{r) 



(5.3) 



which imphes, on using again the definition (2.9), that 0(r) obeys the second-order equa- 
tion 



(f l{l + l) 



dr'' 



+ 



+ P^Eir) 



0(r) 



2 ^4^1 






-0(r) 



(5.4) 



where we have defined (cf Eq. (5.48) in [16]) 



P^^]Ei^)^PwA^)'rll 



w 



" + W'^ (// + (Ai - W^)-i) - 2-W' 



(5.5) 



For example, if a potential W of Coulomb type is considered, one finds from (4.11) and 
(5.5) that the limiting form of the eigenvalue equation (5.4) as r ^ is entirely dominated 
by the term proportional to /i. More precisely, in such a limit Eq. (5.4) reduces to 



dr"^ 



^^2 



0(r) = 



which is solved by 



0(r) = r e 



(5.6) 



(5.7) 



An analogous method can be used for 0(r) defined in (2.16), finding a parameter-dependent 
potential 
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(5.8) 



which leads again to the limiting behaviour (5.7) when W{r) = -, but now for 0(r), as 
r ^ 0. We can therefore see, in a physically relevant example, that our approach, leading 
to second-order equations for O and O, recovers qualitative agreement with the analysis 
in [16], where it is shown that, no matter how singular is the central potential at r = 0, 
the Dirac operator is always well defined as long as /i 7^ 0. In other words, our formula 
(5.5) accounts clearly for the dominating effect of the anomalous magnetic moment with all 
potentials diverging at the origin. However, a rigorous result on the relation between our 
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approach and the squared Dirac operators studied in [16] remains an interesting technical 
problem whenever W{r) 7^ (cf [13]). 

6. Concluding remarks 

The contributions of our paper, of technical nature, consist in the application of analytic 
techniques that can help to understand some key qualitative features of central potentials 
for the Dirac equation, with emphasis on the mathematical formulation of relativistic 
eigenvalue problems. Although the methods used in our investigation are well known 
in the literature, the overall picture remains, to our knowledge, original (see comments 
below). In particular, we would like to mention the following points (at the risk of slight 
repetitions) . 

(i) The forms (2.10) and (2.17) of the second-order equations for the radial parts of the 
spinor wave function, with Pw,e{'>^) and Pw,e{'>^) defined in (2.11) and (2.18), respectively, 
is very convenient if one wants to understand whether the potential can affect the self- 
adjointness domain of the free problem. 

(ii) The identification of the domains of (essential) self-adjointness of the operators defined 
in (2.12), (2.19a) and (2.19b) is helpful as a first step towards the problem with non- 
vanishing potential W{r), and clarifies the general framework. 

(iii) A potential of Coulomb type, although quite desirable from a physical point of view, 
leads to some non-trivial features with respect to the non-relativistic case. We have in fact 
seen that Pw,E{r) and Pw,E{r) fail to be infinitesimally form-bounded with respect to the 
squared Dirac operators in the free case, if W{r) contains a Coulomb term. Moreover, 
the limit-point condition at zero for the potential in the second-order operators in the 
interacting case is only fulfilled if the inequalities (4.14) or (4.15) hold. In other words, the 
essential self-adjointness on C^{0,oo) of the second-order operators with non- vanishing 
potential is still obtained, but under more restrictive conditions expressed by (4.14) and 
(4.15). This may have non-trivial physical implications: if essential self-adjointness fails 
to hold, we know from section 3 that different self-adjoint extensions of the second-order 
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operators exist, characterized by the choice of regular boundary condition at r = (cf 
(3.20)). The lowest values of / (for which (4.13) does not hold), corresponding to the 
bound states of greater phenomeno logical interest, might therefore find an appropriate 
mathematical description within the framework of self-adjoint extensions of symmetric 
operators. It remains to be seen how much freedom is left, on physical ground, to specify 
the boundary conditions for the self-adjoint extension. 

(iv) On considering the effect of the anomalous magnetic moment, Eq. (2.10) is replaced 
by Eq. (5.4), with the potential term defined in (5.5). For all potentials diverging at the 
origin, the effect of the anomalous magnetic moment is then dominating as r ^ 0. 

Indeed, as far as the Dirac operator is concerned, one can prove its essential self- 
adjointness on C^fR"^ — {0} j in the presence of a Coulomb potential provided that I7I 

(see (4.11)) is majorized by |v3, as is shown in [16], following work by Weidmann (see 
page 130 in [16] and references therein). In our paper, however, we have focused on 
second-order differential operators, and the consideration of a central potential, with the 
associated Hilbert space 

has eventually led to the second-order operators occurring in (2.10) and (2.17) and acting 
on square-integrable functions on the positive half-line. Our calculations, summarized in 
the points (i)-(iv) above, remain therefore original. We should notice that the condition 
I7I < ^a/3 found in [16] is compatible with our inequalities (4.14) and (4.15) for all / > 2. 
In other words, the condition on 7 ensuring essential self-adjointness of the Dirac operator 
leads also to essential self-adjointness of the second-order operators studied in our paper, 
whereas the converse does not hold (one may find a I7I smaller than \/3 but greater 
than |v3). Our analysis has possibly the merit of having shown that some extra care is 
necessary when / = 0, 1, but this should not be unexpected, if one bears in mind from 
section 3 that already in the free case the values / = 0, 1 make it necessary to perform a 
separate analysis (cf [20]). 
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We should also acknowledge that in [21] the essential self-adjointness of powers of 
the Dirac operator had been proved, but in cases when the potential V is smooth. In 
particular, when the potential is a C°° function on R^, no growth conditions on it are 
necessary to ensure essential self-adjointness of any power of the Dirac operator [9, 21]. 
In our problem, however, we have considered a Coulomb term in the potential, which is 
singular at the origin. Although a regular solution of the eigenvalue problem exists [8, 18], 
since the origin remains a Fuchsian singular point, the domain of essential self-adjointness 
of the second-order operators in the interacting case is changed. This is reflected by the 
inequality (4.13) for the fulfillment of the limit-point condition at zero, which now involves 
7, and hence the atomic number [11, 12]. Note also that, to find a real-valued solution 
which is regular at the origin in a Coulomb field, one only needs the weaker condition 
k'^ > 7^ [II7 12]. Thus, a careful investigation of the essential self-adjointness issue picks 
out a subset of the general set of real-valued regular solutions. 

For simplicity, we have considered in the end of section 4 only one 'linear' term. 
More precisely, however, two linear terms are often studied, of scalar and vector nature, 
respectively [8]. Moreover, a naturally occurring question is whether one can extend our 
qualitative analysis to study the (essential) self-adjointness issue for operators involving 



the square root of the Laplacian [22], i.e. J —c^fi A +mQC'^ — ^^. Such problems have 
been the object of intensive investigations, but more work could be done from the point 
of view of rigorous mathematical foundations. In the light of the above remarks, there is 
some encouraging evidence that new insight into the choice of phenomenological central 
potentials can be gained by applying some powerful analytic techniques along the lines 
described in our paper. In the near future, one might therefore hope to re-interpret from 
a deeper perspective the previous work in the literature, including the class of potentials 
responsible for quark confinement. 
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